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ABSTRACT 



Context. The origin of the disks around Be stars is still not known. Further progress requires a proper parametrization of their structure, 

both spatially and spectrally. This is challenging as the disks are very small. 

Aims. Here we assess whether a novel method is capable of providing these data. 

Methods. We obtained spectro-astrometry around the Pa/3 line of two bright Be stars, a Col and ( Tau, to search for disk signatures. 

The data, with a pixel-to-pixel precision of the centroid position of 0.3. .0.4 milliarcsecond is the most accurate such data to date. 

Artefacts at the 0.85 mas are present in the data, but these are readily identified as they were non-repeatable in our redundant datasets. 

This does illustrate the need of taking multiple data to avoid spurious detections. 

Results. The data are compared with model simulations of the spectro-astrometric signatures due to rotating disks around Be stars. 

The upper limits we find for the disk radii correspond to disk sizes of a few dozen stellar radii if they rotate Keplerian. This is very 

close to observationally measured and theoretically expected disk sizes, and this paper therefore demonstrates that spectro-astrometry, 

of which we present the first such attempt, has the potential to resolve the disks around Be stars. 

Key words, techniques: high angular resolution - techniques: spectroscopic - stars: individual: a Col, f Tau - stars: emission-line,Be 



1. Introduction 

For decades it had been surmised that Be stars are surrounded by 
disk-like structures. At first this was inferred by indirect means 
such as the doubly pea ked Ha emis sion line profiles dStruvd 
1193 lb and polarization (lPoeckerllll97 5). This notion was con- 
firmed only much l ater by direct, interferometr ic observations at 
radio wavelengths JDoughertv & TavlonI 19921) . Later, dedicated 
long-baseline optical and near-infrared (NIR) interferometry re- 
solved the disks at selected baselines (e.g. lOuirrenbach et all 



1 19971 iTvcner et al. "2004": Meilland et al."2007 - for a general re- 
view on Be stars see Porter & Rivinius 2003) . So far, relatively 
few Be stars have been studied in this manner. This is due to 
the fact that the disks are small, even the largest observed disks 
are typically of order a few milliarcsec (mas) in diameter and 
observations remain challenging. 

In this paper we investigate the potential of spectro- 
astrometry to detect disks around Be stars. This technique is 
a powerful tool; it enables us to investigate small scale struc- 
tures with a standard instrumental set-up. In addition, since data 
are taken at high spectral resolution, it also allows kinemati- 
cal studies to be performed at superior resolution than for in- 
terferometry. Spectro-astrometry is a proven method to study 
otherwise unresolved structures in longslit spectra. It has been 
used to study binaries (Bailey 1998; Baines et al 2004, 2006; 
Schnerr et al. 2006), outflows fro m young objects ( Takami et al. 
2003i) . disks around young objects dPontoppidan et al.i20 08) and 



even ma de possible the disco very of bi-polar jets from Brown 
Dwarfs JWhelan et alj l2005h . Conceptually, this technique is 
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straightforward, it measures the relative spatial position of spec- 
tral features from a longslit spectrum. For example, the red- and 
blueshifted emission of a rotating disk will be located on oppo- 
site sides of the continuum. Even when spatially unresolved, the 
centroid position of the spectrum will be offset from the con- 
tinuum, and this can be deter mined very accurately to sub-pixel 
values (see e.g. lBaiIevlll998l) . The method has been shown to 
detect binaries at separations of 0.1 arcsec in conditions where 
the seeing was in excess of 2 arcsec, while brightness differences 
between the binary components of up to 6 magnitudes have been 
observed as well (iBaines et al.l 12006). Observationally, it is a 
comparatively cheap method, requiring only a stable spectro- 
graph and a digital detector. It can therefore be applied to large 
samples of object. 

As for example demonstrated by iTakami et alj (l2003h . the 
positional accuracy of the centroid mainly depends on the num- 
ber of photons and the seeing and can be expressed as o" = 
0.5 X FWHM X N^'-, with the error cr and full width half max- 
imum of the profile (typically the seeing) expressed in arcsec 
or milliarcsec, and A^ is the number of photons. For shot-noise 
dominated statistics, A^"^ equates to the inverse signal-to-noise 
ratio (SNR) of the total spectrum. Therefore, the requirements 
are proper sampling, high SNR, an d a narrow instrume ntal point 
spread function (i.e. good seeing). iBaines et al](l2006l) achieved 
a position accuracy, as measured from the root-mean-square 
(rms) variations in the position, of 2 mas in 2 arcsec seeing. The 
aim of the present study is to significantly improve on this statis- 
tic to investigate whether we can detect the milliarcsecond scale 
disks around Be stars. 

This paper is organized as follows. In Sec. 2 we describe the 
observations of two bright, nearby Be stars, and reduction pro- 
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cedure. In Sec. 3 we present the results of the study, introduce 
the sub-milHarcsecond spectro-astrometry and we discuss the re- 
sults in terms of a simple model. We conclude in Sec. 4. 

2. Observations and Data Reduction 

For this experiment we selected two Be stars that were bright, 
close-by and had a track record of strong hydrogen recombina- 
tion line emission. These factors should ensure that they are sur- 
rounded by comparatively large disks. Indeed, ^ Tau had been 
measured to have an Ha diameter of 7 mas (e.g. iTvcner et alj 
I2OO4I) . a Col has no published interferometric data, estimates 
indicate a lar ger size of its line emitting region than of ( Tau 
dPachs et al.ll 992). For the choice of telescope, we had to trade- 
off between an excellent sampling of the spectro-astrometry and 
the choice of target line. Ha may be expected to form in larger 
regions, but the near-infrared instrumentation described below 
provided excellent sampling. For this pilot study, we decided 
to push the best available equipment to us and observed Pa/? at 
1.28//m. 

The spectro-astrometric data were obta ined in service mode 
with the Phoenix instrument (iHinkle et al.l2 003) mounted on 8m 
Gemini South in Chile during the night of December 19 (UT) 
2004. Phoenix is a high resolution near-infrared spectrometer 
operating in the wavelength region 1-5/im. The target line was 
the Pa/? hydrogen recombination line, a strong line in a region of 
the spectrum that is relatively unaffected by telluric absorption. 

The grating was set such that the 1.28 yum Pa/3 line was in 
the centre. The detector was an Aladdin 1024x1024 InSb array 
with a pixel size of 5.9xl0"^jum (in wavelength, corresponding 
to 1.4 kms"') spectrally, resulting in an unvignetted wavelength 
coverage of 1300 kms '. The pixel size in the spatial direction 
was 85 milliarcsec, and the slit was 4 pixels wide. The seeing 
during the observations was hovering between 0.40-0.55 arcsec 
as measured from the central parts of the spectra. The set-up en- 
sured that both the spatial and spectral resolution elements were 
sampled by 4-5 pixels. This is much better than Nyquist sam- 
pling and having the data sampled by such a large number of pix- 
els is a crucial constraint when dealing with spectro-astrometry. 
Although not strictly necessary with such bright targets at this 
wavelength, the observations were done in the standard a-b-b-a 
nodding on the slit to remove sky emission. 

In addition to the usual East- West and North-South slit po- 
sitions, that makes sure we can measure the position angles for 
any extended material in the data, we observed at the opposite 
angles as well. The rationale is that the observations have to be 
repeatable and spurious effects should be identified by multiple 
observations. Real effects are each other's mirror image (one is 
looking at the objects upside down as it were), instrumental ar- 
tifacts would b e present in th e s ame direction on th e array (see 
for more details ^ailevlll998l and lBaines et al.ll2006l) . In the case 
of a Col, the slit position angles (PA) were set at 0°, 90°, 180°, 
and 270°. ( Tau was observed at 58°, 148°, 238°, and 328°. The 
choice of different PAs rather than the usual EW-NW settings 
was that the sht woul d be aligned with the dis k resolved in in- 
terferometric data of (" Ouirrenbach et al]|1997[) . However, they 
report a position angle of -58°, and the omission of the minus 
sign in our instrumental set-up means that the data are somewhat 
less efficient than they could have been. At each slit position, the 
objects were observed 4 times with the nodding procedure and 
the total spectra consist of 16 integrations with exposure times 
of 6s each for a Col and 10s for ^ Tau respectively. The overhead 
associated with rotating the slit was small and the total on-target 
time was less than 45 minutes in both cases. 



The data were red uced in a sta ndard manner for optical data 
using both the IRAF (lTodvlll993l) and Starlink software pack- 
ages. Dark frames were subtracted from the original frames, 
which were then divided by a normalized flatfield. The inten- 
sity spectra were extracted and the individual spectra co-added 
to arrive at the total spectra that are discussed in the remain- 
der of this paper The observational details and some derived 
parameters are summarised in Table 1 The wavelength calibra- 
tion was performed by identif ying telluric a b sorpti on lines mea- 
sured from the catalogue by iHinkle et al.l (Il995h and finding 
the dispersion of the spectra. The resulting wavelength should 
be accurate to within 0.25 kms"^ The FWHM of telluric lines 
that were assumed to be unresolved was measured to be in the 
range of 6. .6.5 kms ', in agreement with the expected resolution 
/i/Ai- 50,000. 

Spectro-astrometric information was extracted from the 2- 
dimensional longslit data by fitting a Gaussian profile to the stel- 
lar flux at each pixel in the spatial direction. By visually inspect- 
ing the data we confirmed that the Gaussians are a good repre- 
sentation of the point spread function. However, as I Porter et aU 
(2004) pointed out, the precise shape of the fitting function is 
not very critical. The positions of the center were recorded as a 
function of pixel number, and later put on a wavelength scale. 

Four position spectra were obtained at every position angle. 
The individual traces had different shapes on the larger (more 
than tens of pixels) scales. This effect is presumably present be- 
cause they were recorded at different locations on the array be- 
cause of the nodding procedure. To eliminate these large scale 
fluctuations, while at the same time preserving the smaller scale 
properties, the position spectra were fitted by a high order poly- 
nomial. The position data taken at opposite position angles (0- 
180° etc), were subtracted from each other to minimize any re- 
maining instrumental effects. Prior to combining, all eight traces 
per orientation were visually inspected in order to identify and 
remove outlying data. For both orientations of a Col, one trace 
each was discarded, for ^ Tau all traces at a PA of 58° were 
retained, while 2 of the 8 traces perpendicular to this were ex- 
cluded. 

The precision in determining the photo-centre of the longslit 
spectra as measured from pixel-to-pixel deviations from the 
mean position was of order 0.3 mas (see Table 1). This positional 
accuracy is the best quality data ever published from a statisti- 
cal pixel-to-pixel point of view. However, in the current data set 
we also see variations, unrelated to the pixel-to-pixel variations, 
over larger scales, with excursions up to 1 mas. The largest such 
multi-pixel variations can be seen in the East- West direction for 
a Col (see Figure [TJ. A sine-wave type feature spanning many 
pixels with an amplitude of 0.85 mas, almost 3cr, appears visi- 
ble. The dip corresponds to the blue peak of the Pa^ line profile, 
but the two local "maxima" around the dip do not correspond 
to any obvious feature in the emission line. Further inspection 
of the data revealed that the individual traces taken at opposite 
angles also show this pattern. Significantly however, the minima 
and maxima occur at different wavelengths, and are therefore not 
reproduceable. A similar behaviour is also found in the data of 
^ Tau, but in this case it is cancelled out after combination of 
the individual traces. We therefore conclude that the traces in a 
Col that seem to display a "periodic" signal with an amplitude 
of less than a hundredth of a pixel are artefacts, as they are not 
reproduceable in our, redundant, data. The reasons for this are 
unclear and future observations are planned to investigate this 
issue. In summary, the data have a, statistical, precision of 0.3 
mas, while larger scale artefacts of order slightly less than 1 mas 
are identified. 
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Object 


HD/HR 


Sp type 


V 


vsin/ 


PA 


exp time 


SNR 


Pa/3 


EW 


Av 


rms pos 








(mag) 


(kms"') 


C) 


(s) 




(A) 




(kms"') 


(mas) 


a Col 


37795 
1956 


B71Ve 


2.64 


176 


90-270 
180-360 


8x6 
8x6 


1200 


-8 




134 


0.35 
0.40 


f Tau 


37202 
1910 


B4IIIpe 


3 


310 


58-238 
148-328 


8x 10 
8x 10 


1500 


-5 




226 


0.25 
0.30 



Table 1. The targets a nd observational details. The fundamental parameters are taken from the Bright Star Catalog 
dHoffleit & Jaschek '1982'). A v denotes the peak separation of the line profiles. The SNR is measured in a line free region of the total, 
co-added, spectra. The rms in the positions are also measured in line free regions, this time in the co-added data in each orientation. 



3. Results 

3.1. The Pap lines 

The results for a Col and ^ Tau are plotted in Fig. [T] The top 
panels display the total intensity spectra while the bottom two 
panels represent the spectro-astrometry at the two orientations 
respectively. Let us first discuss the Pa/J profiles. Both stars have 
doubly peaked emission lines, a Col has a regular, symmetric 
line profile, while the blue peak of ^ Tau's emission is much 
stronger than the red peak. The Equivalent Widths {Wa) are -8 
and -5 A and the line peak separations are 134 and 226 kms"' 
for a Col and f Tau respectively. 

The emission at the wavelengths covered by the Pa/3 line is 
due to three components. Firstly, the emission line itself, sec- 
ondly continuum free-free emission and thirdly, the stellar con- 
tinuum, which is diluted by the underlying, photospheric, Pa^ 
absorption line. We can estimate the contribution of the free- 
free emission to the total flux. At shorter wavelengths it is 
fairly low, as the continuum excess due to free-free emission in- 
cre ases towards longer wa velengths (se e e.g. the in-depth stud 
by lDoughertvetal.lll99ll) . For ( Tau, 'Dough erty eTall (11991 



derive an excess of 0.1 magnitude at 1.25 /im, which roughly 
corresponds to 1 0% o f the emission being due to the disk. 
Dougherty et al. (| 1 99 lb did not include a Col in their sample. 
Dachs et all (1 19881) observed the object in the optical and near- 
infrared one month apart and we derive a quasi-simultaneous 
V-J colour of -0.18. According to Koornneef (1983), (V-J)q 
for a B7V object is -0.25 mag (he does not list values for sub- 
giants with luminosity class IV). Taken at face value, we would 
thus obtain a non-physical, negative excess. However, the dif- 
ference is comparable to the uncertainty in spectral class and 
photometric errorbars, and illustrates that the excess continuum 
emission at Pa/? must be very small. The depth of the underly- 
in g photospheric Pa/? absorption can be assessed using the data 
of [Wallace et al.l (12000). They present medium resolution spec- 
troscopy of 88 MK spectral standards, amongst which a num- 
ber of B-type stars. The stars in their sample with spectral types 
closest to ours, B7III, B7V and B3IV (HR 1791, HR 3982 and 
HR 6588 respectively), span a wide range in spectral type. The 
central dip ranges from 0.6. .0.65 of the continuum for the nar- 
rower fines (B3IV, B7III) to 0.78 (B7V) for the broader line. At 
-100 and -f 100 kms ' from the line center, where the line emis- 
sion peaks are found, the absorption for all three objects reaches 
down to 0.85. .0.9 of the continuum, i.e. a depression of 10.. 15%. 

Hence, given that the spectral types of our target objects of 
B4III and B7IV, are similar to that sampled by these MK stan- 
dard stars, we assume that the photospheric absorption line un- 
derneath the line peaks is about the same fraction of the line free 
continuum. As the peak line emission is roughly twice that of 
the stellar continuum (1.98 and 1.84 for a Col and ( Tau respec- 



tively), we find that the emission from the stellar photosphere 
and the hydrogen recombination line are approximately equal. 

3.2. The spectro-astrometry 

The spectro-astrometric traces, at both orthogonal orientations, 
are shown in the middle and bottom panels of Fig. [T] The traces 
are normafized to the stellar continuum and the deviations from 
it are expressed in milliarcseconds. Several things are immedi- 
ately apparent from the data. Firstly, the rms variations around 
the mean position are much smaller than 0.5 mas (Table 1). Thus, 
the data have a sub-milliarcsecond precision, and therefore con- 
stitute the most accurate spectro-astrometry of any object hith- 
erto observed. As we are exploring unknown territory, it may 
not come as a surprise that we encounter new problems in the 
data. After the multi-pixel variations that stretch to slightly less 
than 1 mas in a Col, the second obvious finding in the spectro- 
astrometric traces is that the telluric features show strong, in fact 
the strongest, signals. These narrow absorption lines are unre- 
solved and we suspect that pixellation effects give rise to these 
features. Checks on data taken at opposite angles revealed that 
the amplitude of the excursion varies with the location of the 
longslit spectrum on the array. It is always in the same direction, 
and therefore an artefact. For illustration we kept the data in the 
figures. Finally, there is no obvious excursion in position space 
associated with the Pa^S profiles. 

We conclude that in the present data we find no evidence 
for significant features in the positional data down to sub- 
milliarcsecond levels. It also illustrates the need for multiple ex- 
posures to ensure consistency and to avoid artifacts in the data 
being interpreted as real. In the following we will investigate the 
implications for the presence of disks around both objects and 
the presence of a binary companions. 

3.3. On the stars' binarity and their disks' sizes 

Based on radial me asurements 4 " Tau is reported to be a single- 
lined, close binary ('Ja radlfl^STh . The separation is 5 mas and, 
from the mass function, th e primary is at lea st 5 magnitudes 
brighter than the secondary (Tyc ner et al.ll2004l) . This large mag- 
nitude difference combined with the small sepa ration makes 
detection of the secondary very difficult (see also iBaines et al.l 
2006), and explains why we do not see a binary signature in the 
data of ^ Tau. a Col has not been reported to be a binary, and the 
data do not show the presence of a binary companion either. 

We can estimate the size of the Pa/3 line emitting region from 
the spectrum. Most methods such as the W^ of the emission em- 
ploy the entire line profile to do this (e.g. iGrundstrom & GiesI 
2006), even interferometrically determined sizes are based on 
the total line emission. Here, we exploit the fact that we have 
spatial information available at high spectral resolution. If the 
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Fig. 1. Data on a Col (right) and ^ Tau (right). The top panels represent the total intensity spectra. The bottom panels are the 
photo-centers expressed in mas, at perpendicular orientations whose values are indicated in the figures. 



disks are in Keplerian rotation, we can compute the distance 
from the star of the bulk of the orbiting material using its rotation 
velocity. Using va lues for the masses and rad ii for the spectral 
types (taken from IStraizvs & Kuriliend [19811 and interpolated 
between B3 and B5 to arrive at a value for ^ Tau), we computed 
the Keplerian rotation speeds at the stellar surface (457 and 489 
kms ' for a Col and ^ Tau). The observed velocities at the line 
peaks (half the peak separation in Table 1) combined with the 
distances to the objects provided by Hipparcos, yield the dis- 
tance of the line peak forming regions from the starof8.8mas(a' 
Col) and 3.9 mas (^ Tau). However, the observed velocities are 
smaller than the true value by a factor sin/, and the distance from 
which the emission originates is smaller by (sin/)^ (for Keplerian 



rotation). Taking the inclinations derived bv lFremat et al.l ( I2005L 
45° and 66°) we obtain 4.4 and 3.2 mas for a Col and ^ Tau re- 
spectively. The amplitude of the excursion in the positional data 
associated with these separations can be calculated by simulat- 
ing the data at the line-peak, convolving them wit h the seeing 
and determining the spectro-astrometric trace (cf. IBaines et alj 
2006). In both cases, the photo-center is precisely halfway be- 
cause the line peaks are equally bright as the star. We therefore 
would expect, based on the above, that in the present data line 
excursions of up to 2.2 mas {a Col) and 1.6 mas (^ Tau) can 
be observed. This is the maximum observable separation, as the 
above computation assumes that all line emission arises from a 
thin ring with a rotation speed corresponding to the line peak. 
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Fig. 2. Toy model predictions of the spectro-astrometric signa- 
ture of a disk surrounding a Be star with properties similar to 
the target objects. The top panel shows the flux spectra, the bot- 
tom panel shows the predicted astrometry for the set-up used in 
our observations for three different cases, an outer radius of 70 
stellar radii (solid line), 30 stellar radii (short-dashed), 10 stellar 
radii (long-dashed). 



Not all emission at the line peak comes from a single ring 
however, as the projected velocities for smaller, faster, rings 
will be observed as well at the observed Doppler shifts. To as- 
sess this eff'ect we performed some simple model calculations. 
We assume the star to be suiTounded by a geometrically thin, 
Keplerian rotating disk reaching onto the stellar surface, with 
the line flux per unit area following a simple power law in ra- 
dius. The main input parameters of the model are the stellar ra- 
dius, rotation speed, the inclination and emission line strength 
(which are all fairly well known), the remaining free parameters 
are the disk's outer radius and the exponent of the power law. 
The model produces a two dimensional position-velocity dia- 
gram, which is binned up and smoothed to represent our pixel 
sizes of 5 kms"', 85 mas and seeing of 500 mas respectively. 
From the resulting data, the spectro-astrometry is then mea- 
sured. Changing the outer radius of the model disk increases 
the spectro-astrometric excursions, which then occur at lower 
velocities, as expected from Keplerian rotation. A stronger line 
flux will yield a larger spectro-astrometric excursion because the 
photo-centre shifts more in the direction of the emission line. 

In the extreme case of optically thick emission, the power- 
law will have a flat slope and the emerging line flux will be 
dominated by the outer parts of the disk. In the other extreme, 
that of optically thin emission, the power law depends on the 
density distribution. For an isothermal, flaring Keplerian disk, 
the surface density, and by implicatio n the flux per unit area, 
has an r"^ powerlaw dependence (cf. ICarciofi et al.ll200 o). As 
a consequence, the line emission moves towards the inner parts 
of the disk. The main positional excursions will thus occur at 
higher velocities, closer to the star and therefore be smaller than 



in the optically thick case. Changing the exponent of the pow- 
erlaw also affects the shape of the emission line. A shallower 
exponent, more representative of the optically thick case, puts 
more flux at lower velocities, while a steeper power law, closer 
to the optically thin situation, results in naiTower lines, with the 
line peak at higher velocities. 

In general though, unless the exponent gets too steep, the 
excursions are of similar magnitude when the same line-to- 
continuum ratio is simulated. We performed a large parameter 
study, but for the purposes of this paper, we will restrain our- 
selves to one illustrative example representative of both objects. 
We set the line-to-continuum ratio to be 2 (as per the spectra in 
Fig. [Hand derived above), use a stellar rotational velocity of 475 
kms ' (halfway the values for both objects) and an inclination 
of 55° (also roughly halfway the two objects) and use a stellar 
radius of 0.2 mas. For the outer radius of the disk we take 10, 3 
and a maximum of 70 stellar radii (cf. Marlboro ugh et al.]|1997h . 
The resulting data are shown in Fig. |2] The top panel presents 
the resulting model line profiles. As expected, the lines are dou- 
bly peaked with peak separations that are larger for smaller disk 
radii. The separations range from ~1 lOkms"' for the largest disk 
to 160 and 270 kms"' for the smallest disks, respectively. This 
trend is explained by the fact that these velocities coiTespond 
to the Keplerian rotation speeds at the maximum possible radii, 
where most of the line flux originates if the emission is opti- 
cally thick. The most notable differences between the model line 
profiles and the observed line profiles are the relative narrow- 
ness of the line peaks and the little emission at low projected 
velocities. This is probably due to the fact that the model disks 
are assumed to be geometrically thin, resulting in low projected 
emitting surface areas at low velocities. In reality the disks are 
flared, and therefore the emitting area will be much larger, in par- 
ticular at these low velocities. In addition, line broadening is not 
taken into account here. For a proper treatment, radiative transfer 
models such as those bv lCarciofi & Bjorkmanl (12006) will be an 
excellent tool. Using such advanced models is beyond the scope 
of this paper, in which we wish to obtain a rough figure for the 
excursions only. 

We also note that the blue peak of ^ Tau is much stronger 
than the red peak. This is most likely due to one-armed oscilla- 
tions in its disk, which give rise to such asy mmetry, as for exam - 
ple detected for ^ Tau at the 0.7 mas level bv lVakili et"an (Il998l) . 
Accordingly, the positional excursion of the blue peak would be 
larger than that of the red one, but its signature would not affect 
the overall appearance of the spectro-astrometry. 

Moving to the spectro-astrometric signature in the simu- 
lations, we find that the positional excursions are smaller for 
smaller disks (1.5, 0.6 and 0.2 mas respectively) for the same 
reason as that the peak separations are larger: most flux comes 
from the outer parts, and the largest model disk will naturally re- 
sult in the largest detection. As the orientations of the disks are 
not aligned with the slit positions, we might in reality observe 
smaller excursions in our data by up to a factor of 0.7 due to 
projection effects. 

The bottomUne of the exploratory model simulations is that 
disks with a size of order 70 stellar radii and a predicted excur- 
sion of 1.5 mas could just about have been observed at the 3-4 cr 
level, while the 30 stellar radii disk would have been a 2 cr de- 
tection. Thes e limits are ap proaching the real sizes of the disks. 
According to lTvcner et al.l (12004). the disk of ^ Tau has a diam- 
eter of 7 mas, and thus a radius of ~18 stellar radii, whereas a 
Col's dis k is assumed to be larger, mainly because of its larger 
line EW dPachs et alJll992l) . It is clear that our pixel-to-pixel 
precision, of order 0.35 mas, is reaching that needed to detect 
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Keplerian disks. It will thus be possible to measure the disks 
and their kinematics and therefore to properly constrain the disks 
with future data. In order to achieve this, observations should re- 
sult in a better SNR, or be targetted at stronger emission lines 
(in terms of line-to-continuum ratio), either for stars with larger 
disks or from intrinsically stronger lines at different wavelengths 
such as Ha. 

High precision data such as these combined with the latest 
radiative transfer models lCarciofi & Biorkmanl(l2006l) will allow 
us to be in a position to fully constrain the kinematical structure 
of Be star disks, and reveal their origin. 

4. Conclusion 

In conclusion, we employed high precision spectro-astrometry 
to assess the potential of the method to detect the disks around 
two Be stars. We achieved rms variations in the position spectra 
of order 0.3 mas, the highest precision spectro-astrometric data 
in the literature. We did not detect any features related to the Pa/3 
lines, but found artefacts at the 1 mas level. These were easily 
identified as they were non-repeatable in our redundant datasets. 
Simple, robust, estimates of the size of the line forming re- 
gions showed that the current set-up was on the limit of detect- 
ing the disks, if they were rotating Keplerian. Indeed, the method 
has the potential to distinguish between Keplerian rotating disks 
and angular momentum conserving disks which would be much 
smaller. This study, the first of its sort, has shown that the method 
has great potential in probing small, sub-milliarcsecond, scale 
structures and future observations using improved set-up, even 
higher SNR, and possibly moving to the intrinsically brightest 
hydrogen recombination lines which should reveal even larger 
disk-signatures in the data are planned. 
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